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Abstract. In the present work a two-component dark matter model is studied adopting the 
(— I ' degenerate scenario in the R-parity conserving NMSSM. The gravitino LSP and the neutralino 

^ NLSP are extremely degenerate in mass, avoiding the BBN bounds and obtaining a high 

I ■ reheating temperature for thermal leptogenesis. In this model both gravitino (absolutely 

stable) and neutralino (quasi-stable) contribute to dark matter, and direct detection searches 
for neutralino are discussed. Points that survive all the constraints correspond to a singlino-like 
neutralino. 
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^ 1. Introduction 

There is accumulated evidence both from astrophysics and cosmology that about 1/4 of the 
I energy budget of the universe consists of so called dark matter, namely a component which is 

■ non-relativistic and neither feels the electromagnetic nor the strong interaction. For a review 

on dark matter see e.g. [T]. Although the list of possible dark matter candidates is long (for a 
nice list see [2]), it is fair to say that the most popular dark matter candidate is the lightest 
super symmetric particle (LSP) in supersymmetric models with R-parity conservation [3j. For 
^ I supersymmetry and supergravity see [3]. The simplest supersymmetric extension of the standard 

^ ■ model that solves the mu problem [5] is the next-to-minimal supersymmetric standard model 

- - - (NMSSM) [6j. If we do not consider the axion [7] and the axino [8], the superpartners that have 

the right properties for playing the role of cold dark matter in the universe are the gravitino 
and the lightest neutralino. By far the most discussed case in the literature is the case of the 
neutralino (see the classic review [9]), probably because of the prospects of detection. However, 
in the case is which neutralino is assumed to be the only dark matter component, one has to 
face the fine-tuning problem and the gravitino problem [10]. In most of the parameter space 
the neutralino relic density turns out to be either too small or too large [11]. Furthermore, 
unstable gravitinos will undergo late-time cascade decays to a neutralino LSP. These decays will 
destroy the light element abundances built up in BBN, unless Tr < 10^ GeV [T2], which poses 
serious difficulties to the thermal leptogenesis scenario [13] . If, on the other hand, gravitino is 
the LSP and therefore stable, playing the role of cold dark matter in the universe, it is then 
the neutralino that will undergo late time decays into gravitino and hadrons, and the gravitino 
problem is re-introduced [Hj . 

It has been shown that in the degenerate scenario [TB] the BBN and CMB constraints are 



avoided, and high values of the reheating temperature are obtained compatible with thermal 
leptogenesis. Here we focus on the scenario in which the masses of the gravitino LSP and 
ncutralino NLSP arc extremely degenerate in mass. Under this assumption neutralino becomes 
quasi-stable participating to the cold dark matter of the universe together with gravitino, and 
it is still around and it can be seen in direct detection searches experiments. 

This article is organized as follows. In the next section we present the theoretical framework. 
In section 3 we discuss all the relevant constraints from colliders and from cosmology, and we 
show our results. Finally, we conclude. 

2. Theoretical Framework 

2.1. The NMSSM 

The particle physics model is defined by the superpotential 

W = eij (Yu Hi Q'u + Yd Hi d + YeH{ U e) - e^^A S HlH^ + ^kS^ (1) 
as well as the soft breaking masses and couplings 

- Aoft = m| Q* Q + m-^ u*u + m?^ J* d + m | L* L + m| e* e 

+m%^ H* Hi + mjj^ H*H2 + m| S*S 
+eij (a„ Yu Hi Q'u + Ad Yd HI d + Y, H\ LP e + H.c.) 
+ [-e,jXA^SH\Hl + \kA^ + H.c.) 
-i (M3 A3 As + M2 A2 A2 + Ml Ai Ai + H.c.) (2) 

When the singlet acquires a vaccum expectation value, S, we obtain an effective fi parameter, 
fieff = A5. Imposing universality at the GUT scale, a small controllable number of free 
parameters remains, namely 

tanP = Vu/vd,mo, Ao,mi/2, X, Ak 
and the sign of the effective /v, parameter can be chosen at will. 

Because of the extra singlet superfield, in the NMSSM there is a larger higgs sector and a 
larger neutralino sector. The neutralino mass matrix is characterized by the appearence of a fifth 
neutralino state, meaning that the composition of the lightest neutralino has an extra singlino 
contribution 

X? = NnB^ + NuWi! + NisH^I + N^H^ + iVis^ (3) 

In the following, neutralinos with Nf^ > 0.9, or Nf^ > 0.9, will be referred to as bino- or 
singlino-like, respectively. 

Furthermore, in the Higgs sector we have now two CP-odd neutral, and three CP-even neutral 
Higgses. We make the assumption that there is no CP-violation in the Higgs sector, and therefore 
the CP-even and CP-odd states do not mix. We are not interested in the CP-odd states, while 
the CP-even Higgs interaction and physical eigenstates are related by the transformation 

h'a = SabH'b (4) 

where S is the unitary matrix that diagonalises the CP-even symmetric mass matrix, a,b = 
1,2,3, and the physical eigenstates are ordered as irif^o < m^o < m^o. 



2. 2. Production of gravitinos 

In the usual case (not in the degenerate scenario) gravitinos can be produced after inflation 
in two ways. One way to produce gravitinos is with scatterings from the thermal bath, and 
another is from the out-of-equillibrium decays of the NLSP, which decouple from the thermal 
bath before primordial Big-Bang Nucleosynthesis and decay after the BBN time. Thus, imposing 
the WMAP bounds [16] we can write for the gravitino abundance 

0.1097 < J^3/2/i^ = f^3/^/i^ + yi/2^^h^ < 0.1165 (5) 

where 

' mNLSP 

with rriMLSP the mass of the NLSP, and ^NLSph'^ the abundance the NLSP would have, had it 
not decayed into the gravitino. The thermal contribution is given by (approximately for a light 
gravitino, <^ rrig) [IT] 
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In the limit where mf^LSP ~^ itt-g and tj^lsp ^ 10^^ sec the scenario looks as if one would have a 



two-component dark matter with the NLSP contribution Qtvlsp^^j and a gravitino contribution 
from thermal production only, . Therefore, in the degenerate scenario with uinlsp — 
the WMAP bound becomes 

(8) 



0.1097 < n^dmh"^ = ^NLSph^ + ^lf2h^ < 0.1165 
where from now on the NLSP is the lightest neutralino, x = NLSP. 



3. Constraints and results 

- Spectrum and collider constraints: We have used the computer software NMSSMTools [18], we 
have performed a random scan over the whole parameter space (with fixed /i > motivated by 
the muon anomalous magnetic moment), and we have selected only those points that satisfy all 
the theoretical requirements, as well as the LEP bounds on the Higgs mass, collider bounds on 
SUSY particle masses, and experimental data from B-physics [HI [20]. For all these good points 
the lightest neutralino is either a bino or a singlino, and contrary to the case where neutralino 
is the dark matter particle, here we do not require that the neutralino relic density falls within 
the allowed WMAP range. 

- As we have already mentioned, the total dark matter abundance, and not the neutralino 
one, should satisfy the cold dark matter constraint [16j 

0.1097 < ncdmh"^ = ^x^^ + ^3/^^^ < 0.1165 (9) 

that relates the reheating temperature after inflation to the gravitino mass as follows 

Q.n = A{m^,my)TR + Q.^h^ (10) 

For a given point in the cNMSSM parameter space, the complete spectrum and couplings have 
been computed, and we are left with two more free parameters, namely the gravitino mass 
and the reheating temperature after inflation. The gravitino mass is equal essentially to the 
neutralino mass, and the precise value can be determined if we specify the neutralino lifetime. 
In the discussion to follow we have used a neutralino lifetime r = 10^^ sec, although the results 



are not sensitive to it, and the figures we have produced for different values of the hfetime cannot 
be distinguished. Finally, the reheating temperature after inflation is obtained from the cold 
dark matter constraint. The thermal production contribution cannot be larger than the total 
dark matter abundance, and for this we can already obtain an upper bound on the reheating 
temperature 

Q / \ /TeV " 

Tr < 4.1 X 10^ 



rrin 



GeV 



(11) 

100 GeV, it 



100 GeV. 

Assuming a gluino mass rrig ~ 1 TeV, we can see that for a heavy gravitino, 
is possible to obtain a reheating temperature large enough for thermal leptogenesis. 

- For neutralino NLSP in the degenerate scenario, the only decay mode is x ^ iG, for which 
the decay width can be computed once the supergravity Largrangian is known [21j, and it is 
given by [a [22] 
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where M* is the Planck mass 
limit where the mass difference Am 
Am <^ 171^, niQ, 



171^ is the neutralino mass, and 9^ 



the neutralino lifetime becomes 



is the weak angle. In the 
is much lower than the masses themselves. 



1.78 X 10^3 



|A''ii cos 9w + -^^12 singly I 
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(13) 



From this formula one can see that for a mostly bino-neutralino a mass difference of 1 MeV is 
already enough to give a neutralino lifetime larger than the age of the universe. 

- Neutralino-Nucleon spin-independent cross-section: LHC is now running and collecting data. 
Although LHC is a powerful machine to look for physics beyond the standard model, it is known 
that other facilities are also needed to offer complementary information towards the direction 
of searching for supersymmetry and identifying dark matter. The gravitino interactions are 
suppressed by the Planck mass, and therefore direct production of gravitinos at colliders and/or 
direct detection prospects seem to be hopeless. On the other hand, for a weakly interacting 
neutralino there are existing as well as future experiments that put experimental limits on the 
nucleon-neutralino cross-section. The spin-independent cross-section is given by 



Ami 



vr 



In 



(14) 



where mr is the Nucleon-neutralino reduced mass, m^ = mNm^/ {rriN + ^x)' ^"^^ 

In 



rriN 



2^ ^Tq + 97 •'^G 



q=u,d,s 



q=c,b,t 



(15) 



In the above, f!^Q - 
matrix elements [23] 
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AN) 
=u,d,sJTq ' 



we have taken the following values for the hadronic 



Ap) 

JTu 



0.020 ±0.004. 



Ap) 

JTd 



0.026 ±0.005. 



f!^^ = 0.014 ± 0.003 , f!^} = 0.036 ± 0.008 , 
and Qg is the coupling in the effective Lagrangian 



/^P) = 0.118 ±0.062, 
f^] = 0.118 ±0.062. 



(16) 



^eff = OiiXXQiQi 



(17) 




Figure 1. Neutralino scattering off a nucleon by a neutral Higgs boson exchange. 




Figure 2. Neutralino scattering off a nucleon by a squark exchange. 

where i = 1, 2 denotes up- and down- type quarks, and the Lagrangian is summed over the three 
quark generations. The coupling Uq can be decomposed into two parts, Oq = Uq + aq, where 
the first term is the t-channel exchange of a neutral Higgs (Fig. 1), while the second term is 
the s-channel exchange of a squark (Fig. 2). The expressions for Oq is terms of the masses and 
couplings of the model can be found in [24j . 

Our main results are summarized in the figures below. In Fig. 3 and Fig. 4 we show the 
Nucleon- neutralino spin- independent cross section (in cm^) versus neutralino mass and lightest 
Higgs boson (in GeV) respectively. The blue region corresponds to a bino neutralino, while the 
green region corresponds to a singlino neutralino, and the curves are the current experimental 
limits from CDMS [25j. According to our results the bino scenario is already ruled out, while in 
the singlino case the upper region can be probed by future experiments. In Fig. 5 we show the 
reheating temperature after inflation as a function of the neutralino/gravitino mass. The blue 
region corresponds to a bino, the blue points correspond to singlino, and finally the red points 
correspond to singlino with relatively high values of the cross-section, namely cr^-]\f > 10~^^ cm?. 
The largest values of Tr correspond to a bino, which is ruled out, and for the singlino with 
relatively high values of cross-section we obtain a reheating temperature Tr ~ 5 x 10^ GeV for 
a neutralino/gravitino mass ~ ~ 200 GeV. In the last figure we show the {1710-1711/2) 
plane (mg and mi/2 in GeV) for singlino points with a cross-section larger than 10~^^ cm?, or 
lower than 10"^'' cm^. We see that mo is not larger than 600 GeV, and therefore future direct 
detection experiments cannot probe a region of the parameter space which can neither be probed 
by LHC. 

4. Conclusion 

In the framework of NMSSM, which solves the mu problem, we have assumed that the gravitino 
LSP and the lightest neutralino NLSP are degenerate in mass. Under this assumption the 
neutralino becomes extremely long-lived avoiding the BBN bounds. In this scenario we have a 
two component dark matter made out of the absolutely stable gravitino and the quasi-stable 
neutralino. We have performed a random scan over the whole parameter space keeping the 
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Figure 3. Spin-independent neutralino-nucleon (proton) cross-section versus neutralino mass. 
Shown are the available experimental bounds from CDMS, and the predictions of the theoretical 
model. The blue region corresponds to a bino-like neutralino, while the green points correspond 
to a singlino-like neutralino. 
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Figure 4. Spin-independent neutralino-nucleon (proton) cross-section versus the lightest Higgs 
mass. Shown are the available experimental bounds CDMS, and the predictions of the theoretical 
model. The blue region corresponds to a bino-like neutralino, while the green points correspond 
to a singlino-like neutralino. 

points that satisfy the available collider constraints plus the WMAP bound for dark matter. 
These points correspond to either a bino or a singlino neutralino. We have computed the 
neutralino-nucleon spin-independent cross section as a function of the neutralino mass and the 
lightest Higgs mass, and we find that the bino is ruled out. Then we explored the (m-o — ^^1/2) 
parameter space, and the reheating temperature dependence of the neutralino/gravitino mass for 
the singlino points that correspond to cross section values to be probed by future experiments. 
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Figure 5. Reheating temperature versus neutralino (or gravitino mass). Blue points correspond 
to bino, green points correspond to singlino, and red points correspond to singlino with a cross- 
section larger than 10~^^ cm^. 

600 r . • 




"1-1/2 (GeV) 

Figure 6. The {mQ-mi/2) plane for the singlino points. One color corresponds to a cross-section 
larger than 10^^^ cm^, and the other color corresponds to a cross-section lower than 10~^^ cm^. 
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